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We present a method of computing any one-loop integral in lattice perturbation theory by system- 
atically expanding around its continuum limit. At any order in the expansion in the lattice spacing, 
the result can be written as a sum of continuum loop integrals in analytic regularization and a 
few genuine lattice integrals ("master integrals"). These lattice master integrals are independent of 
external momenta and masses and can be computed numerically. At the one-loop level, there are 
four master integrals in a theory with only bosonic fields, seven in HQET and sixteen in QED or 
QCD with Wilson fermions. 

PACS numbers: 11.15.Ha, 12.38.Bx 
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I. INTRODUCTION 

Lattice regularization provides a non-perturbative for- 
mulation of quantum field theories and allows for the 
approximate numerical evaluation of the associated path 
integral . Nevertheless, perturbative calculations with 
lattice regularization are often necessary, particularly 
when the problem at hand contains several disparate 
scales. In such situations it is sometimes possible to 
separate perturbative short distance physics from non- 
perturbative long distance physics. This separation de- 
pends on the regularization prescription, and in order to 
combine the perturbative results for the short distance 
part with the non-perturbative evaluation of the long 
distance matrix elements, both must be computed within 
the same regularization scheme. A typical example is the 
weak decay of a hadron. The physics associated with the 
high scale Mw can be treated perturbatively by absorb- 
ing it into the Wilson coefficients of a low energy effec- 
tive Lagrangian. The low energy physics is then obtained 
by evaluating the matrix elements of the operators that 
induce the decay. To match the lattice results for the 
matrix elements with the Wilson coefficients, typically 
obtained in dimensional regularization, a perturbative 
calculation in lattice regularization must be performed. 
Other important applications of lattice perturbation the- 
ory include the construction of actions with smaller dis- 
cretization errors |Q , and the extraction of quark masses 
and the strong coupling constant |Q, ^, ^ from lattice 
simulations. 

It is well known that perturbative calculations with 
lattice regularization are difficult. The propagators and 
interaction vertices in lattice regularization are much 
more complicated than their continuum counterparts; 
lattice regularization also badly violates Lorentz invari- 



ance. Consequently, the standard set of tools for one- 
loop calculations in continuum perturbation theory does 
not appear to be useful in lattice calculations. In the 
standard approach to lattice perturbation theory the rel- 
evant loop integrals are therefore evaluated numerically, 
which has several drawbacks: a) the amount of numerical 
computations necessary for realistic calculations is huge; 
b) cancellations between individual diagrams can render 
numerical results unstable; c) the continuum limit, i.e. 
the limit in which the inverse lattice spacing becomes 
much larger than external momenta and masses, has to 
be taken numerically as well. There exist a number of 
techniques to both reduce the amount and enhance the 
precision of the numerical integrations involved. These 
methods rely on Reisz's power counting theorem Q and 
use momentum subtractions to split the loop integrals 
into lattice tadpole contributions and a remainder whose 
continuum limit can be taken naively. The number of lat- 
tice tadpole-integrals is then reduced by exploiting rela- 
tions among them ||, ^, |ll| and techniques for the precise 
numerical evaluation of such integrals have been devel- 
oped 

In this paper we demonstrate that a completely dif- 
ferent approach to lattice perturbation theory is possible 
if advanced methods of continuum perturbation theory 
are applied to lattice calculations. We first show how to 
construct the systematic expansion of any lattice integral 
around its continuum limit by employing the technique 
of asymptotic expansions [E2| developed for continuum 
loop integrals. After the expansion, the original lattice 
integral is expressed as a sum of two different contribu- 
tions: continuum one- loop integrals that are regularized 
by means of analytic regularization and massless lattice 
tadpole-integrals. We then use integration-by-parts iden- 
tities p3 to systematically reduce the number of tadpole- 
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integrals to an absolute minimum. Those remaining are 
called master integrals and are evaluated numerically; 
their number depends on the theory and varies from four 
in any theory with only bosonic fields to sixteen in QED 
or QCD with Wilson fcrmions. 

The advantage of this approach is that the one-loop 
lattice tadpole-integrals relevant for any conceivable cal- 
culation can be expressed in terms of master integrals in a 
process independent way. Therefore, once these relations 
are established and the master integrals are computed 
numerically, perturbation theory on the lattice reduces to 
perturbation theory in the continuum. 

Our paper is organized as follows. After presenting 
our notation, we illustrate the method by applying it to 
massive tadpole-integrals. Although this example is very 
simple, it exhibits all the basic features of the method. 
We then study Feynman diagrams that involve static and 
Wilson fermions. In both cases the loop integrals can be 
expanded around their continuum limit in exactly the 
same way as the bosonic integrals. However, the com- 
paratively more complicated form of the propagator for 
Wilson fermions leads to a larger number of master lat- 
tice integrals. Finally, we use the technique to rederive 
the known results for the gluon self-energy and the static 
and Wilson fermion two-point functions at one loop order 
in perturbation theory. 



II. PRELIMINARIES 

In addition to particle masses and momenta, the result 
of a calculation in lattice perturbation theory depends on 
the lattice spacing a. It is convenient to define dimension- 
less quantities by multiplying all momenta and masses by 
the lattice spacing so that m = amphys, P = opphys, etc. 
We use the standard lattice notation 



f^Y.Pl PM = 2sin^, 
^=1 ^ 



(1) 



throughout the paper. The massive bosonic propagator 
on the lattice is 



1 



(k^ + m?) 



(2) 



Loop integrals contain products of these propagators in- 
tegrated over the Brillouin zone; a typical loop integral 
has the form 



1 



1 



{2-kY (fc2 ^_ ((p::pfc)2 + 



(3) 



where p is the external momentum and d = 4 is the 
space-time dimension. 

Our goal is to construct a procedure for expanding 
lattice integrals around their continuum limit. In terms 
of dimensionless quantities, this limit corresponds to all 



external momenta and masses of the particles becoming 
small. We begin by mapping the integration region in 
Eq.(^ to an infinite volume and define new integration 
variables ry^. 



77^ tan(A:^/2). 



(4) 



In terms of the new variables, the loop integrations in 
Eq.(^ range from — cx3 to ~\-oo. It turns out to be con- 
venient to introduce variables similar to Eq.(^ also for 
external momenta and to rescale m to to = to/2. The 
self-energy integral in Eq.(^) then becomes 



42 / 11 (1 



d 2 



-, -1 



(5) 



where pi = tan(pi/2). 

The representation of lattice loop integrals as in Eq. (||) 
is similar in form to continuum loop integrals, albeit with 
unconventional propagators. The form of these prop- 
agators makes it impossible to apply the basic contin- 
uum techniques for performing loop calculations, such 
as Feynman parameterization and Passarino-Veltman re- 
duction. However, we will show in the next section that 
more advanced techniques, such as recurrence relations 
and asymptotic expansions, are applicable. 



III. EXPANSION AROUND THE CONTINUUM 
LIMIT 

We demonstrate our method by evaluating the massive 
bosonic tadpole-integral. Although this is the simplest 
loop integral, the treatment of integrals with external mo- 
menta and different masses does not pose any additional 
difhculty; the part of these integrals that has nontriv- 
ial dependence on the masses and momenta, can be ob- 
tained by evaluating continuum one loop integrals. Using 
the variables defined above, the massive lattice tadpole- 
integral can be written as 



G{m) 



°° d 



where 



Db{v) = E 



(6) 



(7) 



We are interested in the continuum limit of Eq.(^, 
TO ^ 0. G{m) therefore depends on a small parameter 
and this fact can be used to simplify the calculation of 
the integral in Eq.(^). Unfortunately, a straightforward 
Taylor expansion of the integrand in to is not possible, 
since G(to) is not an analytic function of to. 
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Traditionally, lattice loop integrals with external mo- 
menta are simplified using Rcisz's power counting the- 
orem [0. In this approach, a number of momentum 
subtractions are performed on the integral. These sub- 
tractions split the integral into two parts: a polynomial 
in the external momenta with tadpole-integrals as co- 
efficients, and a remainder that depends non-trivially 
on the external momenta, but whose continuum limit 
can be taken naively. When applied to theories with 
massless fields, this method produces infrared divergent 
tadpole-integrals, and makes it necessary to introduce an 
infrared regulator at intermediate stages of the calcula- 
tion. These infrared divergences illustrate that the loop 
integrals are not analytic functions of the small momenta 
and masses; this is reflected in the appearance of terms 
such as log(TO). 

As is well known from continuum perturbation theory, 
the fact that the Taylor expansion of the Feynman inte- 
gral in a small parameter does not commute with the loop 
integration is not an obstacle for constructing a highly 
practical procedure to perform the expansion before the 
integration; this technique is known as the asymptotic 
expansion of loop integrals |12|. To apply this procedure 
to lattice loop integrals we must introduce an additional 
regulator into the integral in Eq.(||). We use analytic 
regularization, and replace the integral G(m) by 



G(m) 



°° d 

And J IL (l+T^f) L ^wyj 



-1-5 



(8) 



After the regulator is introduced, the expansion of G{fh) 
in TO can be constructed; the limit S can be taken 
at the end. The result for G{rh) is obtained as a sum of 
two contributions: 



G{m) = Gsoft(n^) + Ghard("T-)- 



(9) 



The soft and hard contributions are calculated by apply- 
ing the following procedure to the integrand in Eq.(|^): 

• Soft: assume that all the components of the loop 
momentum 77 are small, rji ^ rh ^ 1. Perform 
the Taylor expansion of the integrand in Eq.(|^) 
in the small quantities rji and m. The expansion 
coefficients in this region are standard continuum 
one loop integrals, regularized analytically. Note 
that no restriction on the integration region is in- 
troduced. 

• Hard: assume that all the components of the loop 
momentum are large, rji ~ 1 3> to and Taylor ex- 
pand the integrand in to. The expansion coeffi- 
cients are given by the functions H({ai};n), de- 
fined in Eq. ( [T3| ) below. These functions are related 
by algebraic and integration-by-parts identities. 

We note that the two regions discussed above do not 
account for all the possible scalings of the components 
of the loop momentum. A potential contribution arises 



from regions where some components of the loop momen- 
tum are soft and the other are hard. Following the above 
logic, the integrand of Eq.(^ should be expanded in the 
small quantities; it is easy to see that in this case the 
analytic regularization as introduced in Eq.(^) does not 
fully regulate the resulting expressions. The contribu- 
tions of these mixed regions are set to zero. To justify 
this prescription we consider another possible choice of 
the regulator. We regulate the loop integrals by including 
a factor (sin^(fci/2))'^ in the measure for each component 
of the loop momenta. The disadvantage of this regula- 
tor is that the resulting soft parts are more difRcult to 
calculate than in analytic regularization. Its advantage 
is that it fully regulates the mixed regions, where some 
loop momentum components are small and others are 
large. The integrals occurring in the mixed region are 
scaleless and therefore vanish. We now introduce both 
the sine-function and the analytic regulator simultane- 
ously to show that the mixed regions do not contribute. 
Since the sine-regulator does not contain any unregulated 
regions, the original integral is recovered by the prescrip- 
tion 



G(M) = lim limG(M,e,(5). 



(10) 



In this prescription the mixed regions do not contribute. 
After showing that for non-zero S both the hard and the 
soft parts are analytic functions of e at e = 0, we can take 
the limit e at the level of the integrand and thus re- 
cover the integrals in analytic regularization. To prove 
analyticity of the hard part, we use the relations between 
the integrals that appear in the hard part to express them 
through convergent master integrals and verify that there 
are no 1/e singularities in the relations of any hard in- 
tegral to the master integrals. We can therefore drop 
the sine-function regulator and ignore the contributions 
of the mixed regions. The analytically regularized in- 
tegral is then recovered and its expansion is constructed 
according to the rules given after Eq. (||) . We now discuss 
how the soft and the hard contributions to the massive 
tadpole-integral are calculated. 



A. Soft part of the integrals 

The soft part of G(to) is easy to obtain. After expand- 
ing the integrand, we arrive at analytically regularized 
continuum one loop integrals. Higher orders in the ex- 
pansion introduce a large number of strongly divergent 
integrals with non-covariant numerators; a consequence 
of broken Lorcntz invariance on the lattice. This is not an 
obstacle, since all such integrals can be easily evaluated 
using the identity 

^ r(a-d/2-E.«.) -i-rT./i 



r(a) 



(11) 
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The calculation of the soft part of integrals that depend 
on external momenta or several distinct masses is equally 
simple: this part is always given by continuum integrals 
in analytic regularization and can be brought to the form 
of Eq.(|T^ using standard methods, such as Feynman pa- 
rameters. 

The result for the soft part of the massive tadpole- 
integral is 



47r^Gsoft 



(m) 



1 



1 — 2 log(m^ 



f ^ + ^ - log(m) ) + 0(m6). (12) 



B. Hard part of the integrals 

We now consider the hard contribution to G(m). Re- 
call that the hard part is obtained by expanding Eq.(||) 
in a Taylor series in to; the resulting expression is 

47r4G(m)i,ard-ff ({1,1,1,1},!) 

-to2 (l + S) iJ({l,l,l,l},2) 

- 4 
771 

+ — {1 + S){2 + S) i/({l,l,l,l},3) + ... . 
The functions iJ({a,j};n) are defined by 



°° d 



— 71 — 6 



(13) 



where and n are integers. In a theory with only bosonic 
propagators, the functions H{{ai}; n) form the full set 
of genuine lattice integrals needed to perform one-loop 
calculations in lattice perturbation theory. We are forced 
to consider H integrals with at ^ 1 for two reasons: they 
appear in the expansion of loop integrals with external 
momenta and/or nontrivial numerators p^ ; in addition, 
they are needed in the reduction process of an arbitrary 
i/-integral to the master integrals. 

We now use integration-by-parts identities to fully ex- 
ploit the algebraic relations between the various integrals 
H{{ai},n). These relations are derived using the fact 
that the integral of a total derivative vanishes in analytic 
regularization: 



= 



1 



L [i + n^^y. 



(14) 



for each value oi ^ — 1, 2, ... d. Computing the derivative 
in Eq. (^4|) and rewriting the resulting expression in terms 
of the functions i/({ai},n), we obtain an algebraic rela- 
tion between integrals with different values of {ai\ and n. 
Another equation can be obtained by partial fractioning, 



i.e. by using the linear dependence of the five "propa- 
gators" in the function iJ({ai},n). The complete set of 
algebraic relations is therefore 

= {n- +Y.^a+ -1)\h, (15) 
= {l + 2a^{a+ - l) + 2{n + 5)n+ a+{a+ -1)\h . 



The conventions are such that the operator increases 
(decreases) the index by one. 

Similar integration-by-parts relations for lattice inte- 
grals were first studied in ||, ^ , where it was shown that 
the entire class of integrals i/({afc};n) can be reduced 
to d master integrals in d dimensions. Here, we neither 
attempt to solve these equations explicitly nor to rewrite 
them in such a form that the reduction of a given in- 
dex is manifest. Instead, we adopt a brute force strategy 
and use computer algebra to explicitly solve the equa- 
tions for a given range of indices. An efficient algorithm 
for solving such recurrence relations has been described 
in |l^. First, a criterion which selects a simpler integral 
out of any two integrals is chosen. Typically, integrals 
with lower values of the indices are considered to be sim- 
pler. The above equations are then solved for a very lim- 
ited range of indices, using Gauss's elimination method. 
The calculation is repeated after supplementing the cho- 
sen set of equations with a few relations involving higher 
index values. By iterating this procedure, the equations 
( p^ can be solved for the entire index range needed in 
a given calculation. The advantage of this brute force 
method is that it immediately generalizes to integrals in- 
volving more complicated propagators (e.g. those of Wil- 
son fermions) or to higher loops. In the continuum, this 
strategy has lead to the solution of problems that would 
have been very difficult to tackle by manipulating re- 
currence relations by hand, for example the four point 
functions of QCD at two loops [p^ . 

Having established that the bosonic case requires four 
master integrals, we can choose those in any way we like. 
It is convenient to choose master integrals that are con- 
vergent in the limit (5 — *■ 0, since their expansion in 5 can 
be obtained by expanding the integrands in Eq.(p3[) in a 
power series in 5. A possible choice that satisfies this cri- 
terion is i?o = i?({l, 1, 1, 1}; 1), = H{{1, 1, 1, 2}; 1), 
ff2 = H{{\, 1, 2, 2}; 1), i?3 = i?({l, 2, 2, 2}; 1). 

It is a fairly easy task to numerically evaluate the fi- 
nite bosonic integrals to relatively high precision. In the 
case of H integrals, this task can be further simplified by 
using the Schwinger representation for the propagators, 
which permits a reduction to onc-dimensional integrals. 
Writing 



1 



1 



(fc2)a r(a) 



(16) 



the integration over the momentum k factorizes and can 
be done analytically, leaving one dimensional integrals 



5 



over products of modified Bessel functions /„(t/2) for nu- 
merical evaluation. In this way, we obtain the following 
one-dimensional integral representations for the master 
integrals i?o-3 defined above: 



continuum limit, including terms up to 0{m^^). We note 
that the expansion converges up to relatively large values 
of the mass, m ^ 2.5, where the Compton wavelength of 
the particle is less than half the lattice spacing. 



dte'^'lo{t/2t-' {Io{t/2)+h{t/2)y 
x[l + {jE+lnt)S] + 0{S^). (17) 



The precise numerical evaluation of these one dimen- 
sional integrals is straightforward. Our results for the 
master integrals are 



Ho = 60.3676829- 16.7691976(5 
Hi = 36.0154101- 1.83118998(5^ 
H2 = 22.1620171 + 4.75278418 S ^ 
H3 = 14.1541117+ 7.32383339(5 



(18) 



Since it is rather simple to obtain the above results, 
the integrals i?o-3 form a convenient basis. However, in 
order to facilitate the comparison of our results with the 
literature, we point out that other choices of the master 
integrals are possible. One possibility is 



^i/({l,l,l,l};0) = ^2^ 
^H({1,1,1,1};1) = 4^2 

^i/({l,l,l,l};2) = -i-2 1n2 



(19) 



1 + I6 7r2 62, 



4i^({i,i,i,i};3)-- ^ 



2(5 



ln2- 



64 TT^ 6, 



The three constants 61.2,3 are closer to traditional choices 
of the basis numbers. For instance, Ref. [|| uses the 
notation Zq = bi, Fq = 167r^62 + 7e and Ylit^ Z\ = 
26 -h 3 TT^ - 384 TT^ 62 3072 tt^ 63. The two constants 5i 
and hi are equal to P\ and Pn of |l^; however, 63 ^ P3. 
We can use the relations between the two sets of basis 
integrals to obtain the constants 61-3. For example, the 
integral -ff (1, 1, 1, 1, 2) is given by 



i7(l,l,l,l,2) 



(3^5 + 5) ^ 
^ 12(1 -f ,5)^'' 
(2(5 + (5^ - 6) 
2(1-H(5)(5 



(3 + 2(5) 
6(1 + (5) 



i?i (20) 



^^ + (rT^^^- 



Using the numerical results for the integrals 77o-3 pre- 
sented above, we then obtain the value of 62 ■ Performing 
similar calculations for the other integrals, we arrive at 



61 = 0.15493339, 

62 = 0.02401318, 

63 = 0.00158857. 



(21) 



Finally, we present the result for the massive lattice 
tadpole-integral by combining the soft and the hard con- 
tributions discussed above. In Fig.|l| we compare the nu- 
merical evaluation of G(m) with the expansion around its 



0.15 




0.05 



FIG. 1: Comparison of the numerical evaluation of the prop- 
agator at the origin G{m) with the expansion around its con- 
tinuum limit. The thick line corresponds to the numerical 
result while the four dashed lines show the expansions in 
m = mphya a to the indicated orders. Recall that m = 2m. 



IV. OTHER PROPAGATORS 

The scalar lattice propagator in Eq.(||) arises in the 
simplest of infinitely many possible discretizations of the 
continuum action. Unfortunately, this simplest choice 
leads to large discretization errors. These can be reduced 
by introducing improved actions which include higher di- 
mensional operators chosen to cancel lattice artifacts to a 
given order in the expansion around the continuum limit 
0. In doing so, the free propagator usually gets modi- 
fied. A famous example is provided by fermions on the 
lattice, where the simplest discretization leads to the well 
known doubling problem. Among the many solutions 
to this problem, the most common is to follow Wilson's 
suggestion and add the dimension-five term ipOip to the 
fermion action [ |Tot . 

Despite the fact that the inclusion of higher dimen- 
sional operators in more refined versions of the lattice 
action changes the propagator, our strategy remains the 
same. The structure of the soft part of the loop integrals 
remains unchanged, since it encodes the behavior at small 
momenta; the hard part, however, becomes a function of 
the couplings of the higher dimensional operators. This 
typically leads to more complicated recurrence relations 
for the hard integrals, but the approach described in the 
previous section is robust enough to solve the resulting re- 
currence relations. We now illustrate this for both static 
and Wilson fermions. 
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A. Static Fermions 

In a bound state of a heavy (m^ 3> ^qcd) and a light 
quark, the heavy quark remains almost static, and the 
probability of finding virtual heavy quark-antiquark pairs 
is small. In this situation the heavy quark can be treated 
as an approximately static color source. The systematic 
expansion around the static limit is obtained using the 
Heavy Quark Effective Theory (HQET)||l|l. This effec- 
tive theory provides a method of putting heavy quarks on 
the lattice that does not require excessively small lattice 
spacings. 

The discretized heavy quark propagator is llq] 



1 



(22) 



In the static limit the pole in the propagator occurs at 
A:4 = 0, rendering an ie prescription necessary even when 
working in Euclidean space. Using the variables intro- 
duced in Section II, the static propagator can be written 
as 



1 - »??4 

2rj4 + ie 



(23) 



The soft integrals are obviously the ordinary continuum 
HQET integrals. After regulating the bosonic propaga- 
tor, the hard integrals become 



°? d 

iJ({a,};n,m)= J Uji 



dv^ [Db{v)] 



— n — 6 



(24) 



Again, we wish to use integration-by-parts relations to 
reduce the number of independent hard integrals to an 
absolute minimum. Both the relations obtained by dif- 
ferentiation with respect to 771-3 and the partial fraction- 
ing identity for the bosonic propagator coincide with the 
relations given in the previous section. The integration- 
by-parts identity obtained by differentiating with respect 
to 774 becomes 



{2{n + S) {a+ 
- [(3 + m) m++ + 2aia+] m++} H = . (25) 

The partial fractioning identity for the static propagator 



IS 



{{l-a+) m++ -af} H = 0. (26) 
An additional algebraic identity is obtained by replacing 



=7' ZTT (5(774) . 

rji + ie 774 



(27) 



The principal value contribution to H{{ai}; n, 1) vanishes 
because of the symmetry 774 — > —774, and the contribution 
from the 5- function is independent of 04: 

iJ({ai, 02, as, 04}; n, 1) = H{{ai,a2,a3, 0}; n, 1) . (28) 



Performing the loop integration over 774, it is easy to 
see that the integrals H({ai, 02, as, 0}; n, 1) reduce to the 
bosonic hard integrals in three dimensions. The lattice 
HQET master integrals are therefore given by the four 
bosonic master integrals discussed in the previous sec- 
tion, and the three bosonic master integrals in d = 3 . 
For the three-dimensional bosonic integrals we choose the 
basis 

1 



TT 
1 

An 
1 

16 TT' 



^H({1,1,1},0) = 1 + 0((5), 



3 H({1, 1, 1}, 1) = Si + 0(6) w 0.252731010 , (29) 
■ H{{1, 1, 1}, 2) = .S2 + 0(6) w 0.012164159 . 



B. Wilson fermions 

Wilson fermions [|oj are often used to introduce 
fermions on the lattice, and we therefore consider this 
case in detail. The propagator of a Wilson fermion is 



G(fc) 



Eu 7m sin(fc^ 



r 1.2 
2 



(30) 



For our discussion only the form of its denominator, 

4i?F(fc,r7i) =^sin2(fc^) + (77i + ^fc2)2^ (3^) 

is relevant. In the continuum limit, external momenta 
and masses tend to zero while the Wilson parameter 
r 7^ stays fixed. In the notation of Section II the de- 
nominator is 

2 

+{m + rDB{il)f . (32) 



Drill, m) = J2]i 



We consider a general one loop integral with both bosonic 
and fermionic propagators and regulate the fermionic 
propagator. The integrals that appear in the hard part 
of the expansion are 

(33) 



H = H{{ai},n,m) = 



0) 



-m — 5 



We again begin by writing down the algebraic and 
integration-by-parts identities for the hard integrals. We 
find three identities by partial fractioning: 



(34) 



l + ^(a+-l)n+|i/ = 0, 

^ i 

-(r^-l)[5:(a+-l) 

i 

i-+^a+(a+-l 



m 



H ^0, 



H = 0- 
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d additional ones arise from integration- by-parts: 

1 1+2 (m+S) a+(a+-l)m+ 2 a+-l~2 '^{a+-l 
^ j 

+ 2 Gi (a+ - 1) + 2 n a+ (a+ - 1) n+ jiJ = . (35) 

These equations are considerably more complicated than 
the ones encountered in the scalar case. Not only does 
the index space become six-dimensional, but the recur- 
rence relations also shift indices by up to three units. 
While the solution of the bosonic recursion relations in 
the index range needed for the applications of Section 
was accomplished within hours, about a day of com- 
puter time was required to perform the reduction in the 
fermion case for r = 1. We are left with sixteen mas- 
ter integrals, which must be computed numerically. For 
purely fermionic integrals this number reduces to ten. 

Because of the complicated structure of the fermion de- 
nominator, we cannot use the Schwinger representation 
( |l6| ) to arrive at one-dimensional representations for the 
master integrals. A general strategy to obtain accurate 
results for more complicated cases is to choose a basis of 
integrals that are very convergent in the infrared, such 
as those with negative powers of propagators. Because of 
the simple behavior of the integrand at small values of k, 
these integrals can be obtained by replacing the integra- 
tion with a sum over a small lattice. We illustrate this 
point for the purely bosonic integrals. Instead of the four 
integrals Hi defined in Eq. (^^, we could have chosen to 
evaluate four integrals from the set 



n+S 



ff({l, 1,1,1},^ 



(36) 



to second order in S. Replacing ki — > ^{h ~ ^) and 
trading the integral over the Brillouin zone for a sum 
over li — 1, . . . , L, we obtain the value of the above in- 
tegral up to finite size effects of order 1/L^+^"'. Even 
the choice n = 0,1,2,3 gives rather accurate values for 
lattice sizes L > 50. In reference this strategy was 
applied to fermionic tadpole-integrals and precise values 
for the corresponding master integrals were derived. We 
have verified the first few digits of their numerical results. 



V. APPLICATIONS 

To test this method in a practical application, we use 
it to rederive the known results for the gluon and fermion 
self-energies in QCD. The calculation of these quantities 
is instructive since it exhibits all the complications that 
arise in lattice perturbation theory. For example, the ex- 
pansion of the gauge action around the continuum limit 
leads to infinitely many vertices. When they multiply di- 
vergent integrals, sub-leading vertices contribute in the 
continuum limit and lattice calculations consequently in- 
volve a larger number of graphs, than the corresponding 



calculations in the continuum theory. Furthermore, the 
vertices in lattice gauge theory have a complicated form. 
It takes roughly one page to write down the Feynman 
rule for the four gluon vertex The tensor structure 
of the graphs can therefore be quite lengthy. Summing 
over the internal indices in the first graph in Fig. one 
ends up with several hundred terms. However, the large 
number of terms is not a significant obstacle since the 
technique described here can be easily automated. 



(a) 



(b) 





(e) 



(d) 



^(r(r(r6yrij\Mi)jj)jJL 
(f) 



(g) 

FIG. 2: Feynman diagrams contributing to the gluon self en- 
ergy at 0{g^). The graphs d, f and g have no continuum 
analog. The graph g arises from the expansion of the Haar 
measure. 



A. Gluon self-energy 

The individual gluon self-energy graphs of Fig. |^ di- 
verge quadratically with the lattice spacing, and we have 
therefore expanded them to sub-leading order. They con- 
tain terms proportional to the the lattice invariants g^jy, 

Pi, and g^v . In the sum the quadratic divergences 
and the non-covariant pieces cancel and the result for the 
inverse gluon propagator becomes 



+ Fsoft) 



9 



r(rh, 



ard 



(47r)2 

i^+l)SabP^P.. (37) 
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The gluonic contribution to the self-energy is 



C. Static quark self-energy 



p(G) I 5 14 5 „2 



, , 1 5 1 p2 



4ard-^<i 35 9 ^" ■ Ar2 9 



14 



2 7r^ 14 TT^ 



6i- 



^ln4 + e(-^4-^6,+8.^6.-ln2 



where iV = 3 is the number of colors. The infrared di- 
vergence at i5 = in the hard part cancels the ultraviolet 
divergence of the soft part after the two terms are added. 
For r — 1, the Wilson fermion contribution to the gluon 
vacuum polarization is 

where x = m? /p^ , Lg = 1/S — ln(m^/4), and 



The wave function renormalization of the heavy quark 
is infrared divergent. Regulating the infrared divergence 
with a gluon mass A, the self-energy takes the form 

ShQEt(w) — . f ^„ (Shard + Sgoft) , (39) 

where (x = \^/uj^, Ls = l/S- ln(AV4)) and 



Shard = -Sin Si+UJ < - 



(2-C) 



+1652^2 (2-0 -e- 2 (2-0 ln2L 



(2-e)L5 + 4 + e 



\/l + x + 1 



Upon extracting the mass and wave- function renormal- 
ization from the above expression, we recover the result 
ofRef. ilsll. 



f{x) = Vl + 4a;ln 



VT+4~E + 1 



VI + 4a; - 1 



After combining the hard and soft parts, we reproduce 
the result of Kefs. |l|. 



B. Wilson fermion self-energy 

For r = 1, the self-energy of a Wilson fermion with 
mass m is 



47r)^ V 



(4 7r)2 

The hard part is given by 

4L = -51-435, 



(38) 



- ^-4^ + 13.74 - 2.91 e , 



■^haid ^ 

_ 4 

■^hard ^ 



= + 3.45 - 3.41 ^ ; 



the soft part is 



yiO) _n 



In 



1 + X 



y(2) _ 
^soft ^ 



(4 + {-Ls -l + il + x) In^} 



with X — m? Ip^, Lg ~ 1/(5 — ln(TO^/4). 



VI. CONCLUSIONS 

We have presented a simple technique for calculating 
the expansion of lattice integrals around their continuum 
limit. The method is based on two observations. First, 
using techniques developed for continuum integrals, it is 
possible to systematically expand the lattice integrals in 
a series in the lattice spacing. The expansion is non- 
analytic and requires introducing an intermediate regu- 
larization to make the expansion of the integrands pos- 
sible. After the expansion, the original integral is a sum 
of two distinct contributions. The first one is a set of an- 
alytically regularized continuum integrals which encode 
the small momentum behavior of the original lattice in- 
tegral and depend in a non-trivial way on the internal 
momenta and masses. The remaining part is a polyno- 
mial in the momenta and masses with massless lattice 
tadpole-integrals as coefficients. Within a given theory, 
these tadpole-integrals are process-independent and need 
to be calculated only once. 

We have shown how to efficiently calculate massless 
lattice tadpoles. Integration- by-parts identities relate 
different lattice tadpole-integrals. These relations are 
rather complicated, and their analytic solution is diffi- 
cult. However, it is not necessary. A robust approach is 
to use the algorithm suggested in ||l^ which allows the 
solution of the recurrence-relations and the reduction of 
all the lattice tadpole-integrals to a few master integrals 
to be performed in an automated fashion. We have ap- 
plied this algorithm to tadpole-integrals in gluodynamics, 
HQET on the lattice and QCD with Wilson fermions. We 



9 



have illustrated the flexibility of the method by comput- 
ing a number of two-point functions in lattice perturba- 
tion theory such as the gluon and fermion self-energies 
and the static fermion self-energy in an arbitrary gauge. 

In the future, it will be interesting to see if this method 
can be generalized to higher orders in perturbation the- 
ory. It is clear that a similar separation of the integrals 
into process-dependent soft and universal hard parts is 
possible; it should also be possible to derive intcgration- 
by-parts relations for the hard parts of the lattice inte- 
grals. A potential problem could arise from the size and 
the complexity of the system of equations governing the 
reduction of the massless tadpole-integrals at the two- 
loop level. In principle, the size of the system is not very 
relevant, but for practical reasons, such as available CPU 
time and memory requirements, it might impose severe 



constraints. 

Although it may sound paradoxical, more precise data 
from lattice simulations will not eliminate the need for 
perturbation theory. In many cases the nonperturbative 
results from simulations need to be matched by equally 
precise perturbative calculations in lattice regularization 
in order to become phenomenologically relevant. Using 
the method discussed in this paper, such calculations are 
not much more difficult than perturbative calculations in 
continuum field theory. Hopefully, this simplification will 
persist at higher orders. 
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